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NTWO - A NITROGEN PROPERTIES PACKAGE 
by K. J. Pew, R. C. Hendricks, and R. J. Simoneau 
Lewis Research Center 

National Aeronautics and Space Administration 
Cleveland, Ohio 

ABSTRACT 

A family of FORTRAN IV subroutines was developed to calculate the 
thermodynamic and transport properties of molecular nitrogen from the 
triple point to 300° K at pressures ranging from 0. 1 to 200 atmospheres. 
The master subroutine NTWO calculates the remaining state variable 
when given any two of pressure, density, and temperature as input. In 
addition, any or all of the thermodynamic and transport properties - 
enthalpy, entropy, specific heats (Cp, Cv), sonic velocity, viscosity, 
thermal conductivity, surface tension - can be obtained by the state vari- 
able call or by a separate call. A call to NTWO with pressure and en- 
thalpy will also generate all the other properties. 

The NTWO package was developed to be used in conjunction with heat 
transfer and fluid flow calculations. Because of this intent, the package 
was written as a family of independent subroutines with a separate sub- 
routine for each specific property or set of related properties. Thus a 
user with limited needs can work with only those subroutines required 
for his program. 

NTWO - A NITROGEN PROPERTIES PACKAGE 

Cryogenic nitrogen is inert, inexpensive, and generally available in 
large quantities at many research centers. It is used for cooling equip- 
ment, heating or cooling other fluids, as a modeling fluid, and in many 
cases as the primary test fluid in heat transfer and fluid dynamics re- 
search. Many of the authors who have submitted research papers to the 
Cryogenic Engineering Conference have used cryogenic nitrogen as their 
working fluid. 
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NTWO is a FORTRAN IV family of subroutines which was developed 
on a 7094-7044 DCS machine to determine the thermodynamic and trans- 
port properties of nitrogen and be used in conjunction with other research 
programs. The NTWO property package is subdivided into subroutines 
and functions for computing the individual properties. The individual 
subroutines may be used independently of subroutine NTWO if the user’s 
requirements are of a more restricted nature (say, memory or time 
limited). Appendix A is a descriptive breakdown of all the functions and 
subroutines used by the master subroutine NTWO, which is listed in ap- 
pendix B. 

The basis for the program is the work of Strobridge (ref. 1). 
Strobridge curve fit the available nitrogen data with a modified 3eneaict- 
Webb-Rubin (BWR) (ref. 2) equation of state. This work represented a 
major advance in establishing an equation of state for cryogens. The ex- 
cellent hydrogen data of Roder and Goodwin (ref. 3) were fit using a form 
similar to that of reference 1. In the early days of work with cryogenic 
hydrogen such an equation represented a milestone to those working in 
research and development. 

The NTWO property package differs from that of Strobridge (ref. 1) 
in that it was developed to be used in fluid flow and heat transfer calcu- 
lations. As such, there are independent calls for the three state vari- 
ables pressure, density, or temperature (see OPERATIONS SHEET). 

In addition, temperature and all the other properties can be obtained as 
a function of pressure and enthalpy, which is of considerable value in 
forced convection studies. 

While enthalpy and entropy are available in reference 1, the specific 
heats (Cp, Cv), sonic velocity, viscosity, thermal conductivity, and sur- 
face tension were not computed. One of the major reasons that the spe- 
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cific heats were not given is the significant deviations in d ?/3T at near 
critical temperatures which generally lead to very large values of Cv, 
see figure 1. After much effort in this area, we finally decided to use 
numerical techniques to determine aU/aT = Cv and aH/aT = Cp. 

The viscosity and thermal conductivity computational method were 
adapted from Svehla (ref. 4) for the gaseous regime and from Thodos 
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(refs. 5 and 6) for the near-critical regime. The surface tension calcu- 
lation was adapted from reference 7. 

Sample property calculations are found in figures 2 to 13. The plots 
of Cp and Cv exhibit some irregularities in the near- critical regime. 
Smoothing techniques could be applied, however, for most applications 
these are unwarranted because the values returned are within acceptable 
tolerance. The properties P, p, T, H, S are in good agreement with 
those tabulated by Strobridge (ref. 1). The values of k, p, and a are 
limited by the state of the art. For example, the spikes in thermal con- 
ductivity that are known to exist in the critical region cannot presently 
be computed. See reference 8 for discussion of this. 

At the time of writing, two papers (refs, 9 and 10) were brought to 
the attention of the authors. Pierce (ref. 9) prepared a FORTRAN 17 
program for nitrogen and helium for the pressure range of 1 to 2500 
psia (0. 068 to 170 atm), temperatures of 120° to 4500° R (66. 5° to 
2500° K) for nitrogen and 35° to 4500° R (19. 4° to 2500° K), for helium. 
A superficial examination indicates that the essential differences are, 
the flexibility and versatility of NTWO is much greate. , and in refer- 
ence 9 the calculation of specific heats uses d P/d T throughout and 
the surface tension is not calculated. The tables presented in refer- 
ence 9 have not been checked as of this time. Reference 10 represents 
a theoretical work with tabulated results for temperatures above 300° K. 
As 300° K nitrogen is considered noncryogenic, the results should not 
alter the results presented herein or those of reference 1. 

OPERATIONS SHEET 

The nitrogen PVT surface has been subdivided into three regimes: 
two-phase or saturation (KR=1); liquid (KR=2); and gas and/or fluid 
(XR=3), see the sketch on table I. In the two phase regime, KR=1, only 
the saturated liquid and saturated vapor properties are returned through 
COMMON /N2/. In order to avoid the choice of which properties should 
be used to compute quality, its calculation is left to the user. 
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Specification of KS 

The values of KS represent what basic thermodynamic input has 
been supplied to NTWO and what basic thermodynamic information is 
desired as output (see appendix B). 

KS=1 implies that temperature and pressure are specified and 
density is required. NTWO computes density 

P = f (T, P) 

KS-2 implies that temperature and density are specified and NTWO 
is requested to compute pressure 

P = f(T, p ) 

KS=3 implies that pressure and density are specified and NTWO 
is requested to compute the temperature 

T = f(P, p) 

KS=4 is a particularly useful call for heat transfer and fluid dy- 
namic calculations. Given the enthalpy and the pressure, NTWO re- 
turns temperature and density 

T = f(P, H) 

P = f(P, H) 

Specification of KP 

To obtain thermodynamic properties other than P, T, or p, and 
transport properties, various values are assigned to KP (see appen- 
dix B). However, a value must still be assigned to KS, to inform NTWO 
of what state variables are known. 

KP=i. in combination with a value of KS informs NTWO that en- 
thalpy must be determined. 

KP=2 , with KS specified, requests NTWO to compute entropy. 

KP=4, specify KS, requests NTWO to compute specific heat at con- 
stant volume and constant pressure, the ratio of specific heats, and the 
sonic velocity. 

KP=8, specify KS, NTWO returns a value for dynamic viscosity. 
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K?=16, specify K3, NTWO returns a value for thermal conductivity. 

K?=32, specify KS, NTWO returns a value for surface tension. 

To obtain any combination of the above properties, simply add the 
XP values. Thus if all the properties are required, specify KP=63. If 
everything except surface tension is required, set X?=31. If enthalpy, 
entropy, and viscosity are required, specify KP=11, etc. 

KU Specification 

KU=2 is a units specification which gives units as listed in table I. 
KU=1 will return those units of table I, except density will be returned 
as gram-moles /liter. KU=1 represents the working units of the pro- 
gram (see units of ref. 1). Additional room is available to specify other 
conversion factors; consult the authors on this matter. 

The subroutine NTWO uses the call vector and COMMON to com- 
municate with the subroutines within the NTWO package and the users pro- 
gram. The call vector and the COMMON are explained on the operations 
sheet, table I. A sample problem is provided in table I to illustrate the 
calling procedure. Table II is provided to illustrate assembly of the pro- 
grams for machine operation. 

Experience has taught us to check KR when working near the satu- 
ration loci and in the near critical region. Our calculations sometimes 
put us in the two-phase regime when we do not want to be; NTWO, how- 
ever, thinks we are asking for saturation values and sets KR=1. Since 
specifying KR=1 is a code instructing NTWO to calculate saturation 
values directly without determining the regime, saturation values could 
be returned for all subsequent calls. Hence, the user must rer’et KR/1 
for nonsaturation values. 

The NTWO properties package is currently available through the 
authors and will become available through the NASA-software program. 


SYMBOL LIST 


Fortran symbol* 


c 

c 

sonic velocity, cm/sec 


CL 

sonic velocity, saturated liquid, cm/sec 


CVP 

sonic velocity, saturated gas, cm/sec 

Cp 

CP 

specific heat at constant pressure, J/(g)(X) 


CPL 

saturated liquid specific heat, Cp, J/(g)(X) 


CPV 

saturated vapor specific heat, Cp, J /(g) (X) 

Cv 

CV 

specific heat at constant volume, J/(g)(X) 


CVL 

saturated liquid specific heat, Cv, J/(g)(X) 


CVV 

saturated vapor specific neat, Cv, J/'(g)(K) 

H 

H 

enthalpy, J/g 


HV 

enthalpy, saturated vapor, J/g 


HL 

enthalpy, saturated liquid, J/g 

k 

K 

thermal conductivity, tV/(cm)(X) 


KL 

thermal conductivity, saturated liquid, «V/(cm)(X) 


KV 

thermal conductivity, saturated vapcr, *V/(cm)(K) 


XP 

thermodynamic and transport properties specifi- 
cation 


KR 

thermodynamic region specification 


XS 

state relation specification 


KU 

units specification 

P 

P 

pressure, atm 

S 

S 

entropy, J/(g)(K) 


Symbols used in the individual subroutines are identified in tne work 
statement of that subroutine (see appendix A). 
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SV 

entropy, saturated vapor, J/(g)(K) 


SL 

entropy, saturated liquid, J/(g)(K) 

T 

T 

temperature, X 

y 

GAMMA 

ratio of specific heats, Cp/Cv 


GAMMAL 

ratio of specific heats, saturated liquid 


GAMMAV 

ratio of specific heats, saturated vapor 

a 

CSGMA 

surface tension, dyne/cm 

M 

CMU 

dynamic viscosity, g/(cm)(sec) 

P 

D 

Q 

density, g/cm 


DV 

, 3 

density, saturated vapor, g/cm 


DL 

9 

density, saturated liquid, g/cm 


i “r -v » *•' * 


n 



N 
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APPENDED A 

FUNCTIONS AND SUBROUTINES CONTAINED IN SUBROUTINE NT/70 

If tne user wishes to disassemble NTWO into its family of subroutines, 
reference 11 should be consulted. 

BLOCK DATA 

1. Equation of State Coefficients 

2. Saturated Vapor Laws Coefficients 

3. Specific Heat at Zero Pressure Coefficients 

4. Conversion constants for use with other units 

SUBROUTINE 

DCHECK(KU, D) 

This subroutine converts the density D to gram- moles /l iter an: 
checkes to see if D is out of range. The units are specified by KU. 

If D is out of range, the program 

(a) tags the word if Tagged Arithmetic is used 

(b) writes an out of range note and continues 

PCHECK(KU, KR, P) 

Converts the pressure P to atmospneres and cnecks for out of range. 
If KR=1, P is checked for out of saturation range. Units are speci- 
fied by KU. 

Out of range cause P to be tagged (if tagged aritnmetic is used) or 
effects an out of range write out. 

TCHECK(KU, KR, T) 

Converts the temperature T to degrees Kelvin and cnecks for out of 
range. If KR=1, T is checked for out of saturation range. 

FUNCTION SOLVE(XI, F, DF) 

Newton-Raphson iteration given an initial estimate XI, tne function F, 
and its derivative DF. If the convergence has not been attained in 
100 iterations, the result is tagged (or a note is outputed). 
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SUBROUTINE ROOT(XO, X2, FOFX, FUNC, XI) 

Solve for XI, such that FUNC(Xl)=FOFX, where XI lies between 
XO and X2. This is a modified half- interval search tecnnique and 
permits only one root between XO and X2. If a root has not been 
found in 100 iterations, the result is tagged (or a note is outputed). 

SUBROUTINE ROOTX - same as SUBROUTINE ROOT - must be in- 
cluded to avoid recursion 

SUBROUTINE 5PLINA(X, Y, NX, T, NT, YINT, KFD, KERROR) 

A spline curve fit routine used for interpolation 
SUBROUTINE DENS(KU, T, P, D, DL, DV, XR) 

Compute the density D, given the temperature T and pressure P. 
The units are specified by XU. If KR is either specified as 1 or 
returned as 1, the saturated liquid and vapor densities, DL and 
DV, respectively, are computed as a function of T only. 

SUBROUTINE PRESS(KU, T, D, P, KR) 

Compute the pressure P, given temperature T and density D. 
Units are specified by KU. If KR is eitner specified as 1 or re- 
turned as 1, P will be the saturation pressure computed as a func- 
tion of T only. 

SUBROUTINE TEMP(KU, P, D, T, KR) 

Compute the temperature T, given pressure P and density D. 
Units are specified by KU, If KR is either specified as 1 or re- 
turned as 1, T will be the saturation value computed as a function 
of P only. 

FUNCTION D3F(DS1 

Function used to solve for density DS given temperature and pres- 
sure 

DSF = P equation of state - P given 
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FUNCTION DDSF(DS) 

Derivative of the above function and used to solve for density D3 
given temperature and pressure 

DDSF = 9 ^ DSF ^ 

SUBROUTINE PSSS(PSS) 

Compute the saturation pressure PSS given temperature 
FUNCTION TSS(PS) 

Function used to compute saturation temperature given the pres- 
sure PS 

FUNCTION TSSF(TSS) 

Function used to solve for saturation temperature TSS given 
pressure PS 

TSSF = Vapor Pressure Equation - log^ 0 (?3) 

FUNCTION DTSSF(TSS) 

Derivative of the function used to solve for saturation temperature 
TSS given the pressure PS 

DTSSF = - 8 iX S .§. F ) 

3T 

FUNCTION TSF(TS) 

Function used to solve for the temperature TS given pressure PS 
and density 

TSF = Equation of State - PS 
FUNCTION DTSF(TS) 

Derivative of the function used to solve for temperature TS given 
pressure PS and density 

axsF . MISS) 


3T 


SUBROUTINE TEMPH(KU, P, H, T, D, DL, DV, KR) 

For a given pressure P and enthalpy H, this routine will return 
values of temperature T and density D. If XR is eitner set 
equal to 1 or return as 1, saturated liquid and vapor densities DL 
and DV, respectively, will be returned. 

FUNCTION T3HF{TS) 

A function used in conjunction with TEMPH. This routine calls 
DENS with the given P and trial values of T to obtain values of 
density D. The enthalpy routine ENTH is then called and a value 
of H is computed. 

SUBROUTINE ENTH(KU, KR, T, ?, D, H, HL, HV) 

Compute enthalpy H given the temperature T, pressure ?, and 
density D. The units are specified by KU. The region is speci- 
fied by KR. If KR is specified as 1, then the saturated liquid and 
vapor enthalpies HL and HV, respectively, are computed as a 
function of T only. 

SUBROUTINE HSLV(PS) 

Compute the saturated liquid and vapor entnalpies given pressure PS 
and temperature 

HV from HSS, and HL=HV-T*HS3LV 
FUNCTION HSS(PS, DS) 

Compute the enthalpy in the region No. 3 or saturated vapor enthalpy 
given pressure PS, density DS, and temperature 


H = H° + 


Cp 0 dT + (Z - 1)RT + 


0 



FUNCTION HINTF(TS) 



Function used to compute the integral in the enthalpy computation 
based on temperature TS 


T 


12 


\ 


FUNCTION HDINT(DS, DSL) 

Compute the integral used in the enthalpy computation from density 
DSL to density DS. 



FUNCTION HDINTF(DS) 


Function used in the computation of enthalpy 


2 

P 




FUNCTION HSSLVF(PS) 


Function used to compute saturated liquid enthalpy from saturated, 
vapor enthalpy or saturated liquid entropy from saturated vapor en- 
tropy given pressure PS and temperature 

HSSLV = — AV 
3T 

SUBROUTINE ENT(KU, KR, T, P, D, S, SL, SV) 

Compute entropy S given temperature T, pressure ?, and den- 
sity D. Units are specified by KU. The region is specified by 
KR. If KR is specified as 1, the saturated liquid and vapor en- 
tropies SL and SV, respectively, are computed as a function of 
T only. 

SUBROUTINE SSLV(PS) 

Compute the saturated liquid and vapor entropies given pressure 
PS and temperature 


SV from SSS, and SL=SV- HSSLV 
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FUNCTION 333(PS, DS) 


Compute the entropy in region 3 or the saturated vapor entropy 
given the pressure PS, the density DS, and temperature. 

t r p r 

S=S° + f Cp 0 d(ln T) - R ln(HTp) + / - - —(—) dp 

J T i |_P p 2 V aT /T 


FUNCTION STINTF(TS) 

Function used to compute entropy based on the temperature T3 

r TS 

/ Cp Q d(ln T) 


FUNCTION SDINT(D3, DSL) 

Compute the integral in entropy from the density DSL to the den- 
sity DS 


AS = 


l d p 


p 2 v%j T 


FUNCTION SDINTF(DS) 

The integrand of the function used to compute entropy 


1 / 3P 


p 2 

SUBROUTINE SPCHTV(KU, KR, T, P, D, CV, CVL, CVV) 

Compute the specific heat at constant volume CV given tempera- 
ture T, pressure P, and density D. ' The units are specified by 
KU and the region by KR. If KR is specified as 1, the saturated 
liquid and vapor specific heats, CVL and CVV, respectively, are 
computed. 

cv AU . *0 - (P/Pl) 

AT AT 
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SUBROUTINE CVPS(KVP, XR, CVS) 

Compute the specific heat given temperature pressure and density. 
The region is specified by XR. This routine is used to determine 
the incremental changes in internal energy U, and enthalpy E for 
computation of specific heats at constant volume and pressure, re- 
spectively. 

SUBROUTINE SPCHTP(KU, XR, T, P, CP, CPL, CPV) 

Compute the specific heat at constant pressure CP, given temper- 
ature T, and pressure P. Units are given by KU and the region 
is specified by XR. If XR=1, the saturated liquid and vapor spe- 
cific heats, CPL and CPV, respectively, are computed. 

SUBROUTINE SONIC (XU, XR, T, D, GAMMA, C) 

Compute sonic velocity C given temperature T, density D, and 
the specific heat ratio GAMMA=CP/CV. Units are specified by 
XU and the region by XR. 

SUBROUTINE VISC(XU, XR, T, D, MU) 

Compute the viscosity MU, given temperature T and density D. 
Units are specified by XU and the region by XR. 

SUBROUTINE THERM(XU, XR, T, D, X') 

Compute the thermal conductivity X, given temperature T and 
density D. The units are specified by XU and the region by XR. 

SUBROUTINE SURFT(XU, XR, T, SIGMA) 

Compute the surface tension SIGMA, given the temperature T. 
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APPENDIX B 
NTWO LISTING 

The following listing of the master program NTWO displays the major logic cnvolving KS and XP. 

crdPcit Tut sr*u HtiAiicits hn.j itu**. oy v,» 

T*ANVM1RI UF P'AtCUUK ^irHV.FN GtvtN T»9PfR4TU«F T* 

P*E5VJRE UrNSIfV U, L.< I91 m 5U‘Y !♦. si ATI LIAT TUNS AH t 
SPCCIHIJ *Y RS. »MIR*WV 4 .»^K k\) I - 4NS**UK V PMfTPtHT|»'S 
AH E y>TUM(ll t>Y HP. IF R‘( IS 'it'j'M ) lU 5P1CIMM AS |, 

PROPERTIES AHt UHPUfi » )l SHI IPAUIA. 

OINFhSllK K^iun. R PC 2021, ftpul l'l .4t»C4( KM 

CUM«M /9 2/RJ,UL,*>2,HL ,rM, >, St , S Y, l V , C VL , LVV . C P, C PL . C P 9 , i.AMM A. 

1L, G4P*AV,C,Cl ,L VP,NU,*UL,4UV , * »Rl » * 9 ,Sl«*4 
A CAL 4L..ILL,miV,A,KL,R9 

DATA UPC 1 /2,i,6, 7. 10 . II ,14,1 4 , i ri , 19 ,22 ,2 ), 2 f* ,/r,lJ,)l, *4, 14,YH, 

134,42,0, 46,47,5 1,51 ,54, 55 ,58,4 4,62 *61/ 

•ATA *rti /4,4, 6, 7,12*13, 14,15*20,21 ,22./),24,?#,n, II* IF- 1T,3«, 
tM.44.A5. 44,47. S2.MtS4.S5 ,60*61 ,62.61/ 

•ATA *PC> /0,4, 10,11. 12 .19.14. 15, 24 ,25. 25 .27, 28, 24, *), 1 1, 40,41,42, 
t4»« 44*45, 44. 55,54,57, 54,59,60,61 +2 ,41/ 

•ATA KPC4 /14, 17, 14,19,20.21. 22, 2), 24.25.26.27, ?•», 24, in, 11 , 

*4 A* 49, 90, 51, 52. 53, 54, 55,56,5 7, 5 4,59, 6Q,*| ,62,6)/ 

«B ID 110,20*30*40 .AS 

C0M4TC OtNSltT 

tf CALL •fft6UU,T,P.0,OL»0V,K«l 4 

M ID M 

IMPUTE MCStUtt 

tC CALI PMSSU U*T*0,P,AJt) 7 

»• to 

CM PUT C flHKRATWt 

tf CALL TMP(KU,P,0,T,W» |( 

■ IN 

CONPVTt TEMP'. NATURE AdO DENSITY GIVEN PNCSSUAE AAO ENTHALPY 

At CALL VtMPMMD,P»MtT,D4»L»»V,HNI I! 

M IP MMM.2II 40,70,40 

COMPUTE ENTHALPY 

MULL lMN(KD*n,T ,P ,»,N.NL ,MV » I? 

!•■■ Mitt 

IP fKMKtim M, 100, AO 

mtmnkm 
m m tit 

CMPiTC CNTUP7 

Mf CALL tftI<«MMUt.P«0,»»lL»S*l 24 

nt m u« i-M 

IPUP-tiClIIII 12C«l30»t2O 
Itt CMtlNMK 
m « i4« 

COMPUTE IPLUPIC MEATS 

130 CALL |PCMTVUU,«.T»P,D.CV,CVL.CWI J4 

CALL SPCHIPtNU.HA ,l,P,CP,CPt,CPtl )# 

IF |ftN,NE*|» GO IU 200 
KAMLA71ATL 
AAHNAV ft P 4/C Vtf 

CALL tCNSl4ll«l,P.t,0L,lt«,ll 4? 

CALL SONtC tNO.M* I.OL .CANMAL «€L I 44 

CALL KMICIAW,M,T.B«.GANNA4.CVPI 47 

■ ■HO 
IN MPUKP/t) 

MM CP, C 9,44444 

CALL KMIC |NU,«N ,1,0, GANNA ,C) *r 

140 M Ml I'UM 

IP EO-NPC1IIII 150,160*1 90 
150 CONTINUE 
•B tO I TO 

CONPUti VI SC04I TV 


16C IP I NR .ME .11 GO TO 21 C 

CALL JtM,t«U»f ,P,’J,'>1 ,JV,1I 
CALL VI1C Ud.KH.l , )l ,*JL> 
CALL VlST.t*U,*« t t,jv,NUVl 
GJ TJ ) TO 

2 IC CALL V14C UU,ftM,T, J,4Ul 
ITO 07 257 I '1,52 

ihhp-hp', «(i 11 m.:,i j 
2* C LHU«C 
GJ K 1» 

COMMIE INHOL (MMCIII 

QM IP MCM.lt MNttt 
Ml NMIMbf^OAiATiftl 


■■ t5t. 


BSSftlt'iS 

PtCT M tit%t.f«t 


, PM Off MOPACC ft MIM 

- tatter 
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M'lWOl MM MWOfKS.Kl' I.i'.D.II.Mf 

' i ' \ \ \ 

i » \ v \ \ 

i ' \ \ \ \ 

1 1 \ ' 

, V \ \ \ 


K !< 0 Unknown, Hum k Kit returned 
K i < I Saturation 


v \ \ Kit 2 Liquid 
\\\\ Kit 3 lias and or fluid 

\\V- Knthalpy. .1 « 

' v v x ^- IXnsity, i; (in'* 

Pressure, atm 
u Temperature, K 


L Thermodynamic and Transport properties 
\ KP 0 Only \\ p, T returned 

i KP 1 H Enthalpy, .! g, (H), (HL), (HVi 

i KP 2 S Entropy, J (g)(Kl, (S),<SLf,(SV) 

• KP 4 Cv Specific heat at constant volume, J/(g)(K), (CV), (CVL), (CVV) 

\ Cp Specific heat at constant pressure, J/(g)(K), (CP), (CPL), (CPV) 

; 7 Ratio of specific heats, Cp/Cv, (GAMMA), (GAMMAL), (GAMMA V) 

\ c Sonic velocity, cm/sec, (C), (CL), (CVP) 

*, KP 8 i, Dynamic viscosity, g/(cm>(sec), (MU), (MUD, (MUV) 

\ KP 16 k Thermal conductivity, W (cm)(K), (K), (KLi, (KV) 

\ KP 32 r Surface tension, dyne /cm, (SIGMA) 

Input specification of independent properties 
KS 1 p -- f(T, P); given T, P find p 

KS 2 P f(T, p); given T, p find P 

KS 3 T f(P,/i); given P, p find T 

KS,4 T ,p f(P, H); given P, H find T.p 

COMMON N2 KU, DL, DV, HL, HV, S, SL, SV, CV, CVL, CVV, CP, CPL, CPV, GAMMA, GAMMAL, GAMMA V, C, CL, 

CVP, MU, MUL, MUV, K. KL, KV, SIGMA 

When working near the saturation locus, check KR/l each tim£ to be assured of nonsaturation calculations. 

At the beginning of the user's program, 

set KU 1 or 2 and 

CALL TAGSCN (tagged arithmetic) 

If TAGSCN is not available on the machine, contact the authors for a modified package. 

SAMPLE PROBLEM 

CALL NTW0(2. 31, 80. . P, 1 254, H, KR) 

3 

The program has been asked to find the pressure P, corresponding to a density of !. 254 g/cm and a 
temperature of 80 K Furthermore, the program will return values for ff, S, Cv, Cp,y, C, g, and k. The 
values of H and P are returned through the call vector and the remaining values are returned through 
the COMMON statement. 


t T r - 120 ti,K/ 

? ^‘^kITn 

§■ /Two-phase 
r y^Vsaturalion* 
“ Liouid 

Kntropy, S 


KP 8 n 
KP 16 k 
KP 32 r 


* For example; KP 1+2 3 will return enthalpy and entropy; KPH + 8H6=25 will return enthalpy, viscosity, 
and thermal conductivity; KP l*2«4 + 8+16 31 returns everything except a; KP=2 will return entropy. 
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TABLE II. - PROGRAM ASSEMBLY 

1. Set up the return of thermodynamic and transport properties 

COMMON /N2/KU, DL, DV, HL, HV, S, SL,SV, CV, CVL, CVV, 

CP, CPL, CPV, GAMMA, GAMMAL, GAMMA V, C, CL, 

CVP, MU, MUL, MUV, K, KL, KV, SIGMA 

2. Set up tagged arithmetic (if not available, see notes in the subroutine 

FORTRAN listing). 

CALL TAGSCN 

3. Establish the units you are working in KU = (specify 1, 2, ). 

4. Add your Program here. 

5. Add NTWO Package, or BLOCK DATA and required subroutines. 

Data 

NTWO Package or BLOCK 
DATA and Subroutines 

Main Program + Subroutines 

KU = 1, 2 

CALL TAGSCN 

COMMON/N2/-etc. 

Program Dimensions, Common, etc. 

Control cards 
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